I've added docstrings, comments and some alterations in logic.

I've also added this catch statement to avoid an issue when the poisson mu is set to zero.

This doesnt feel elegant, and also possibly doesn't match what you'd expect.

I've added docstrings, comments and some alterations in logic.

I've also added this catch statement to avoid an issue when the poisson mu is set to zero.

MULE/packs/tools/fit_functions.py

Line 42 in 255261c

# Collect the position and amplitudes of the finger plot peaks

This doesnt feel elegant, and also possibly doesn't match what you'd expect.

My only concern with this (that is not decoupled from the previous implementation, I believe) is that when you do this, then sigmaq becomes a floating parameter, and that can be problematic when fitting, but I don't have a good solution rn sadly